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Propagation, of short radio waves beyond the horizon is discussed in terms 
of rejection from layers in the atmosphere formed by relatively sharp gra- 
dients of refractive index. The atmosphere is assumed to contain many such 
layers of limited dimensions with random position and orientation. On this 
basis, the dependence of the propagation on path length, antenna size and 
wavelength is obtained. 

INTRODUCTION 

It was pointed out several years ago 1 that power propagated be3'ond 
the radio horizon at very short wavelengths greatly exceeds the power 
calculated for diffraction around the earth. This beyond-the-horizon 
propagation has stimulated numerous experimental and theoretical in- 
vestigations. 2 Booker and Gordon, 3 Yillars and Weisskopf 4 and others 
have developed theories based on scattering of the radio waves by tur- 
bulent regions in the troposphere. This paper proposes a theory in which 
uncorrelated reflections from layers in the troposphere are assumed re- 
sponsible for the power propagated bej'ond the horizon. 

In developing this theory, some arbitrary assumptions of necessity 
have been made concerning the reflecting layers since, at the present 
time, our detailed knowledge of the atmosphere is insufficient. However, 
calculations based on the theory are found to be in good agreement with 
reported measurements of beyond-the-horizon propagation. 

Measurement of the dielectric constant of the atmosphere 5 has shown 
that relatively sharp variations in the gradients of refractive index exist 
in both the horizontal and vertical planes. Although the geometrical 
structure of the boundaries formed by the gradients is not well known, 
one may postulate an atmosphere of many layers of limited extent and 

* This material was presented at the I.R.E. Canadian Convention, Toronto, 
Canada, October 3, 1956. 
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arbitrary aspect.* The number and size of the reflecting layers, as well 
as the magnitude of the discontinuities in the gradient of dielectric con- 
stant which form them, influence the received power. 

The reflecting properties of the layers are discussed first. Next, an ex- 
pression for the received power is obtained by summing the contributions 
of many layers in the volume common to the idealized patterns chosen 
to represent the transmitting and receiving antenna beams, f This ex- 
pression is then used to calculate the effect on received power of changes 
in such parameters as the orientation of the antennas, wavelength, dis- 
tance, and antenna size. 

The MKS system of units is used throughout. 



REFLECTING 
' SURFACE 




Fig. 1 — Reflection by a layer. 



EFFECT OF LAYER SIZE 

Propagation from a transmitting antenna of effective area A T to a 
receiving antenna of effective area A R by means of a reflecting layer is 
illustrated in Fig. 1. The ray from transmitter to receiver grazes the 
layer at angle A. The reflection from the layer depends on the ampli- 
tude reflection coefficient, g, which is a function of the grazing angle, and 
on the dimensions of the layer relative to the dimension of a Fresnel 
zone, t Three cases, depending on the layer dimensions, will be considered. 

* After this paper was submitted for publication, a report was received giving 
some measurements of sharp variations in dielectric constant gradient and esti- 
mates of the horizontal dimensions of layers in the troposphere. J. R. Bauer, 
The Suggested Role of Stratified Elevated Layers in Transhonzon Short-Wave 
Radio Propagation, Technical Report No. 124, Lincoln Laboratory, M.I.T., Sept., 

1956. 

f Under some conditions, layers outside the volume common to the antenna 
beams may contribute appreciably to the received power. Phenomena such as 
multiple reflections and trapping mechanisms are not considered in this study. 

t The power received by reflection from the layer in Fig. 1 can be calculated 
approximately by assuming it to be the same as the power that would be received 
by diffraction through an aperture in an absorbing screen, the dimensions of the 
aperture being the same as the dimensions of the layer projected normally to the 
directions of propagation. The field at the receiver is calculated from the distri- 
bution of Huygens sources in the aperture. The received power, expressed in 
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Case 1 . Large Layers 

If the layer were a plane, perfectly reflecting surface of unlimited ex- 
tent, the power at the terminals of antenna A R would be the same as 
the power received under line-of -sight conditions, 6 

p _ p A T A R 
4X 2 a 2 

If the layer has an amplitude reflection coefficient, q, the received power 

is, 



Pr = Pr 



A T A R '< 



This relation applies when the layer dimensions are large in terms of the 
wavelength an d ar e large compared with the Fresnel zone dimensions; 
that is, b > \/2aX/A and c > \/2aX. 

Case 2. Small Layers 

When the dimensions of the layer are small compared with the Fres- 
nel zone, but large compared with the wavelength, the received power is 
given by the "radar" formula, 

Pr = Pr -r-rr- c (6A) q 

This relation applies when b < \/2a\/A and c < \/2dX. 
terms of Fresnel integrals, is 

Pit = Pt'^TT K7"(«) + S*(v)][C*(v) + SHv)) 

C , bA 

where u = — j= and v = — j=. 
\/\a -y/Xa 

When u and v are very large, we have, approximately, 

C(«) = S(u) = C{v) = S(v) = I 

and the expression for Pr reduces to that given for Case 1 above, except for the 
factor q 2 . 

When both u and v are very small, we have approximately, 

C{u) = a C(v) = v 

S(u) = o S(v) = o 

and the expression for Pr reduces to that given for Case 2. 

When u is large and v is small the expression for Pr given in Case 3 results. 
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Case 3. Layers of Intermediate Size 

If the layer dimensions are such that c is large but 6A is small, com- 
pared with the Fresnel zone dimension, the received power is given by 

2\ s a 6 

In the atmosphere, c and b are likely to be about equal, on the average, 
and we have for this case, \/2aX < b < \/2aX/A. 

All three of these cases may be present at various times, since the 
structure of the atmosphere changes from day to day. However, for the 
purpose of the present study, Case 3 is considered most prevalent and is 
assumed in all the calculations to follow. 

Many of the numerous layers that are assumed to contribute to the 
received power are not necessarily horizontally disposed, they may be 
oriented in any direction. Therefore, reflection in the direction of the 
receiver can take place from layers located both on and off the great 
circle path. If there are N contributing layers per unit volume in the 
region V common to the radiation patterns of the transmitting and re- 
ceiving antennas, then for Case 3, 

2X ,1 a ,J Jv 

In this relation it has been assumed that the layer size and the number 
of layers per unit volume remain sensibly constant throughout the com- 
mon volume. 

The integration process requires expressions for the reflection coeffi- 
cient q and the glazing angle A of the layers in the common volume. These 
quantities are derived in the following sections. 

REFLECTION COEFFICIENT OF A LAYER 

The reflection coefficient of a plane boundary (Fig. 2) separating two 
media whose dielectric constants, relative to free space, differ by an 
increment dc is given by Fresnel's laws of reflection. For both polariza- 
tions, the plane wave reflection coefficient of the boundary is 

q = de/4A- 




c+de 
-i 
Fig. 2 — Reflection at a boundary between two homogeneous media. 
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provided 1 » A" » dr. This reflection coefficient for an incremenfal 
change in dielectric constant can l>e used to calculate the reflection from 
discontinuities in the gradient of the dielectric constant of the atmos- 
phere such as those shown for a stratified medium at y = and \j — h 
in Fig. 3(b). 

Such variations of dielectric constant are assumed to be representative 
of discontinuities in gradient as they exist in the physical atmosphere. 
The variations form the reflecting kiyers. 

The method of calculating the reflection coefficient of such a stratified 
medium is due to S. A. Schelkmioff' and is illustrated schematically in 
Fig. 3 in which the medium has been subdivided into incremental steps. 
Consider the reflected wave from a typical incremental layer, dy, situated 
a distance y above the lower boundary of the layer, 0. From Fig. 3(a) 
it is clear that the phase of this wave is 4.ry/\ sin A relative to that of a 
wave reflected from the lower boundary. The incremental reflection co- 
efficient is </f/4A = —A' (////4A", where A' is the change in gradient of 
the dielectric constant at the boundaries of the layer. The field reflected 
by layer dij is therefore, 



dE R = -E,~e 
4A 2 



-J(47r;//X) 



dy 



One now obtains the complete reflected field by summing the reflec- 
tions from all increments within the layer of thickness h. 



E, = / dE r = jEi 



AX 



H -e 



-j(4WiM) sin Ai 



Jo 16k A' 1 sin A 

This relation shows that the layer is equivalent to two boundaries at 




Fig. 3 — Plunc-wuve reflection al an incremental layer </// within a stratified 
medium extending from y = to // = h. 
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y = and y = h, each with reflection coefficient 

•"IB (2) 

If the abrupt change in slope, the solid line in Fig. 3(b), is replaced 
by a gradual change as indicated by the dotted lines, (2) still holds pro- 
vided d < X/4A. For more gradual changes, d = nX/4A where n > 1, 
the reflection coefficient is 
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Smoothing of the boundaries reduces the value of q. 

It will be assumed in all the calculations to follow that reflection from 
layers in the troposphere is described by (2). 

VARIATION OF STRENGTH OF LAYERS WITH HEIGHT 

The formula for the reflection coefficient includes the factor K, which 
represents the change in the gradient of the dielectric constant at the 
boundaries of the layer. A dielectric constant profile constructed of 
many randomly positioned gradients is shown schematically in Fig. 4. 
The variations are shown as departures from the standard linear gradient. 
Measurements 6 indicate that the fluctuations of the dielectric constant 
normally decrease with height above ground. The changes in the dielec- 
tric constant gradients associated with these fluctuations probably vary 
in a similar manner so that K is some inverse function of the height above 
the earth. However, to simplify the computation of received power, to 
be described later, we have adopted the cylindrical coordinate system 
shown in Fig. 5, and it is convenient, then, to let K be a function of p, 
the distance from the chord joining the transmitter and receiver to the 
point in question. 

We assume, therefore, that 

K = 3 ^22^ (3) 

p 

where Ki is the change in gradient at point A in Fig. 5 which, for a typi- 
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Fig. 4 — Schematic illustration of variation of the dielectric constant in the 
troposphere. 
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Fig. 5 — Coordinate system for a beyond-the-horizon circuit. 

cal path length of 200 miles, is about 1,600 meters (1 mile) above the 
earth or, since p = 211, 3,200 meters from the z axis. 
Equation 3 is used in all the calculations to follow. 

THE GRAZING ANGLE A 

The grazing angle A at the slightly tilted layer shown in Fig. 6 is 
given by 



tan 2A = 



•lap 



fi — 9-2 — n 2 
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Throughout the volume common to the antenna patterns, A « 1, p 
« a and z < a/2. Then 

A « t (4) 

a 

It is evident that A is constant and equal to p/a when the point (p, z) 
is located on a cylinder with axis TR and radius p. It is this feature that 




a 



-a -*| 

Fig. 6 — Grazing angle A at a layer. 




Fig. 7 — Idealized antenna patterns used in this study. 
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dV, = 2(f -a)fl Pi <i fii dv 2 = 2(a-/ + z a )flp 2 dp 2 

Fig. 8 — Integration over the common volume. 
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suggested the unusual idealized antenna patterns shown in Fig. 7, which 
are described in the next section. 

CALCULATION OF THE RECEIVED POWER 

Substituting (2), (3) and (4) in (1), one obtains for the received power, 
P R = P T MaA k A 7 .X"' f p"" dV (5) 



where 



M » 2000b 2 KiN (6) 



To integrate over the volume common to actual antenna patterns 
would be difficult. We have, as mentioned before, replaced the actual 
patterns with the idealized patterns shown in perspective in Fig. 7 and 
in plane projection in Fig. 8. The patterns (Fig. 7) are bounded by side 
planes of the large wedge and by surfaces of cones with axis TR. The com- 
mon volume is indicated by broken lines and is well defined. Since the 
grazing angle A is constant for the incremental cylindrical volumes dVj 
and dV-i shown in Fig. 8, it is easy to integrate over the common volume 
V and we obtain 

£'-"-w'© (7) 

(8) 




The function /(a/0) is plotted in Fig. 9. 

The gain of the idealized antennas is G = 8w/a^(a + 20) and the ef- 
fective area is 

A = aP(a + 20) (9) 

The area of a cross section of the antenna pattern is bounded by two 
straight sides, ra, and two curved sides rdff and r(8 + a)/3. The aspect 
ratio is defined as the ratio of the sum of the lengths of the curved sides 
to the sum of the lengths of the straight sides. It is equal to one when 

B = -£» CO) 

Substituting (10) in (9) gives the effective area of the idealized antenna 
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Fig. 9 — The function / (a/0). 



with aspect ratio one, 



A = i 



(11) 



Substituting (7), (10) and (11) in (5) gives for identical transmitting 
and receiving antennas with aspect ratio one, 



MX 3 1 



X H — * T ~7Z . 



3 a'flWq ,a\d 

" T 1 



/I? 



(12) 



For actual antennas, a may be taken to be the half -power beam-width. 

In the following sections, (12) will be used to derive some of the gen- 
eral properties of propagation beyond the horizon. 

* K. Bullington has suggested that a useful form for equation (12) is 



Pr = 



|_4a 2 a 4 J [_ 30 



'1 



>'(:-)' 



2+r 



where the first term in brackets represents the power that would be received in 
free space, the second term involves the characteristics of the troposphere, and 
the third term is a correction factor for narrow beam antennas. 
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RECEIVED POWER VERSUS ANGLE 

The angle is the angle between the lower edge of the idealized an- 
tenna pattern and the straight line joining the terminals (see Fig. 7). 
The minimum value of is determined by the profile of the transmission 
path. If X, a and a are constants, (12) can be used to calculate the ratio 
of the powers, P Ri and P m , received at two different angles, 0i and 2 , 

,14-- ((- 
P« 2 Wo , a 77^ 



2 +F / [a 



n 



Equation (13) shows the importance of having the angle as small as 
possible. For example, for 6\ = a and 2 /0i = 1-25, the power ratio is 
3.4 (5 db). Thus in an actual circuit the antenna pattern should be close 
to the horizon plane. 



RECEIVED POWER VERSUS WAVELENGTH 



Consider a given path in which a and are specified. Equation (12) 
can be used to calculate the ratio of received powers, P R \ and P R2 , cor- 
responding to two different wavelengths, X] and X 2 . 

Pn, \\J \aj 2j /W 

" J \e 

where m and a-> are the beamwidths of the antenna patterns at wave- 
lengths Xi and X 2 respectively. 

Case I. Equal antenna gains at the two wavelengths. 
For this case, a { = a-> and equation (14) reduces to 



Pr\ /Xi 



Pn2 Vx, 

In free space the power ratio would be 

Phi \\v 

or 

Pr\/Pri (Beyond-Horizon) _ Xi 
PrJPr'i (Free Space) X 2 



(15) 

(16) 
(17) 
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Thus if 400 mcs and 4,000 mcs were propagated over the same path, the 
antenna gains being identical for the two systems, then, on the average, 
one would expect the received power relative to the free space value at 
400 mcs to be 10 db higher than that at 4,000 mcs because of the charac- 
teristics of the troposphere. 

Case II. Equal antenna apertures for the two wavelengths. 
For this case, «i/a 2 = Xi/X 2 and (14) reduces to 

2+? f(i) 

(18) 



Pjn 



Experimental data for this case was obtained on the 150 nautical mile 
test circuit between St. Anthony and Gander in Newfoundland. 8 The 
antennas for both wavelengths were paraboloids 8.5 meters in diameter. 
Simultaneous transmission tests at Xi = 0.074 m and X 2 = 0.6 m were 
conducted for a full year. For this circuit, 

= 0.94° (4/3 earth radius) 

a, = 66 ^ = 0.575° 
8.5 

« 2 = 66 £? = 4.65° 
8.5 

Using these values in (18), we get for the ratio of received powers, 

Pnx/Pm = 1.01 

For antennas of equal aperture in free space, 

Pkl/Pu = (V*i)' = 65.5 

Therefore, 

Pri/Pri (Beyond Horizon) = l_ = _ lg l db 
Pri/Prs (Free Space) 65 

The Summary of Results, Sections 1 and 2 on page 1316 of Reference 8, 
gives - 17 db for this ratio. The agreement between calculated and meas- 
ured values is very good. 
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RECEIVED POWER VERSUS DISTANCE 

If antenna size and wavelength are specified, (f2) gives for two dis- 
tances, ai and 02 , 



(19) 







Phi 
Phi 


9 J. ° 

_ /*y 2 + s 


>d) 








'© 




For a 2 : 


= 2oi. 


■ (19) give 


?s for different values of a/di 




a/0i = 
Prx/Pri - 


= 0.5 

= 276 (24 db) 


1 2 

197 (23 db) 138 


: (21.5 db) 


4 

104 (20 db) 



Fig. 1 in Bullingt on's paper, 9 which gives the median signal level in 
decibels below the free space value as a function of distance, shows an 
18 db increase in attenuation when the distance is doubled. This cor- 
responds to a ratio of received powers of 18 + 6 = 24 db. The examples 
in the table above give an average increase in attenuation of 22 db. 

RECEIVED POWER VERSUS ANTENNA SIZE 

Equation 14 can be used to calculate the effect on received power of 
changing simultaneously the size (and, hence, the beamwidths) of the 
antennas used for transmitting and receiving, the wavelength and dis- 
tance remaining fixed. 



Pn> \di/ „ _i_ £5 r ( a "- 



where P , ( \ and P ff2 are the received powers corresponding to the antenna 
beamwidths a\ and a 2 respectively. 

As an example, let a 2 be constant and equal to 4° and let 6 be 1°, 
corresponding to a 200-mile circuit. The table below gives the ratio 
Phx/Phi as ai is varied. 

a, 4° 2° 1° 0.5° 0.25° 

Pri/Pm (db) 10 18.5 25.7 31.4 

Change in db 10 8.5 7.2 ~" 5/T^ 

Since a is inversely proportional to the antenna dimensions, the table 
shows that continued doubling of the antenna dimensions results in less 
and less increase in output power. The increase varies from 10 to 5.7 db 
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in the tabic. This is a characteristic feature of beyond-the-horizon prop- 
agation. In free space, doubling the antenna dimensions would result 
in a 12-db increase in output power. 

Large antennas and high power transmitters are costly, and a proper 
balance between their costs requires careful studies which are outside 
the scope of this paper. In general, it is not believed worth while from 
power considerations to increase the antenna size much beyond the di- 
mensions that correspond to a pattern angle, a, equal to angle 9. 

Another factor to be considered, however, is the effect of antenna size 
on delay distortion in beyond-the-horizon circuits. From simple path 
length considerations, one concludes that the delay distortion decreases 
when the beamwidths of the antennas are made smaller. Therefore, delay 
distortion requirements may dictate antenna sizes that are not justified 
by power considerations alone. 

SEASONAL DEPENDENCE 

Both the effective earth radius, R e , and the magnitude of the discon- 
tinuities in gradient, Ki , are related to the season of the year. During 
the summer when the water vapor content of the air is high, the effective 
radius and the discontinuities in gradient are larger than in winter. Sub- 
stituting a/R e for 6 and assigning summer and winter values for R e and 
Ki , (12) may be used to calculate the ratio of the power received in 
summer and in winter. 

aR e w r (aRes 
P R (Summer) = / K W Y / R e s\ 5 "*" a J \ a 

P R (Winter) " \K 1W ) \R c w) 2 + «^s ./ aR c w x ^ 

KuY/Rcs 

Ai 17/ \Re W 

For example, if we assume K\ S = 2K XW and R ea =1.2 R c w , then P RS / 
P RW = 11.9 (10.75 db). A seasonal variation has been observed. 8, 10 

DEPENDENCE OF RECEIVED POWER ON ANTENNA ORIENTATION 

The variation of received power with orientation of the antennas at 
the terminals of a beyond-the-horizon circuit differs considerably from 
that observed under line-of-sight conditions. Consider, for example, Fig. 
10 which shows the beams of the transmitting and receiving antennas 
elevated simultaneously. The variation of received power can be calcu- 
lated from (13). As an example, consider the 188-mile circuit between 
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Crawford Hill, N. J., and Round Hill, Mass., for which experimental 
data is published. 10 For this circuit a = 0.05° (3 db points) and di = 1° 
(4/3 earth radius). The table below gives the calculated variation of re- 
ceived power as angle 6-> is varied. 

02 = 1° 1.1° 1.2° 1.4° 1.6° 1.8° 2° 2.2° 
10 log, (P«i/P« 2 ) = 2.3 4.5 8.5 12 15 17.9 20.5 

The received power versus elevation angle, y = 6 2 — 6\ , is plotted in 
Fig. 10. The calculated and experimental curves are in good agreement. 

If the beams of the antennas are steered simultaneously in the hori- 
zontal plane, Fig. 11, the calculation of the variation of received power is 
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Fig. 10 — Relation between received power and vertical angle y. 
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comparatively simple. In horizontal steering, the intersection of the axes 
of the antenna beams moves along line AB in the figure labelled "Cross 
section at 0." If the intersection of the beams moved along the circle 
A-C, the received power would not change. The decrease in power caused 
by moving the beams from position A to B is given by (13). The calcu- 
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Fig. 11 — Relation between received power and azimuth angle 5. 
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lations are identical with the calculations for elevation steering; the ele- 
vation angle 8 2 is related to the azimuth angle 8 and the beamwidth angle 
a by 

8* = (02 - 0i)(0 2 + 0i + «) 

A calculated curve of received power versus azimuth angle is shown in 
Fig. 11 for the Crawford Hill-Round Hill circuit together with the re- 
ported experimental data. 10 The agreement is considered good. 

THE VALUE OF FACTOR M IN EQUATION (l2) 

An average value for the factor M can be obtained from propagation 
data. Using equation (12), the ratio of received powers corresponding to 
free space and beyond-horizon transmission is 

0.75 5 ( 2 +? 



P K (free space) = \ 6 1 ,^ 



P R (beyond-horizon) M ^ , 



i) 



This ratio was found experimentally to be 5 X 10 6 (67 db) for the circuit 
between St. Anthony and Gander in Newfoundland. 8 For this circuit, 
a = 0.081, 9 = 0.0164, X = 0.6. Substituting these values in (22) we 
obtain, 

M = 3 X 10~ 14 (23) 

Substituting this value for M in (12) leads to the following equation for 
a beyond-horizon tropospheric circuit, 

P R = p rX io" 14 - — 2 — — / (j) (24) 

^ + d 

Equations (6) and (23) give 

tfKiN = 1.5 X 10 -17 (25) 

Although values of the layer dimension, b, the change in gradient, Ki , 
and the number of layers per unit volume, N, are not known, it is interest- 
ing to calculate N from (25) assuming reasonable values for b and Ki . 
Assuming Ki = 4 X 10 -8 , which is half the value of K', the average 
gradient of the dielectric constant in the troposphere, and b = 1,000 
(1 km) we find N = 10 -8 or 10 layers per cubic kilometer. 
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CONCLUDING REMARKS 

The interpretation of propagation beyond the horizon in terms of re- 
flection from layers of limited size formed by variations in the gradient 
of the dielectric constant of the atmosphere leads to relatively simple 
results which are in good agreement with reported experimental data. 
The received power depends on the wavelength, the distance, and the 
size of the antennas used for the circuit and on the strength and size of 
the reflecting layers. 

As mentioned earlier, the structure of the atmosphere may change 
markedly from time to time so that large, small and intermediate size 
layers play their parts at different times. Furthermore, the effective size 
of a given layer may be different for widely separated wavelengths, de- 
pending on the roughness of the layer in terms of the wavelength. All 
that can be expected of a study such as the present one is that it serve 
as a guide for estimating the roles of the various parameters involved in 
beyond-the-horizon propagation. 
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